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Bound states for the quantum dipole moment in two dimensions 2 

Abstract. We calculate accurate eigenvalues and eigenfunctions of the Schrodinger 
equation for a two-dimensional quantum dipole. This model proved useful for the study 
of elastic effects of a single edge dislocation. We show that the Rayleigh-Ritz variational 
method with a basis set of Slater-type functions is considerably more efficient than 
the same approach with the basis set of point-spectrum eigenfunctions of the two- 
dimensional hydrogen atom used in earlier calculations. 



1. Introduction 

In a recent paper Dasbiswas et alp] discussed the bound-state spectrum for a straight 
edge dislocation oriented along the z axis. Within the continuum model, the authors 
reduced the problem to the two-dimensional Schrodinger equation for a quantum dipole. 
They obtained very accurate eigenvalues and eigenfunctions by means of a discretization 
of the x — y space. Since this real-space diagonalization method (RSDM) is not 
so practical for highly excited states, those authors also carried out a Rayleigh-Ritz 
(RR) variational calculation with the basis set of eigenfunctions of the two-dimensional 
Coulomb problem [2]. 

The variational RR eigenvalues are known to approach the exact ones from above. 
However, in the present case the ground-state energy e RR = —0.0970 calculated with as 
many as 400 basis functions exhibits a considerable discrepancy with respect to the same 
eigenvalue obtained by RSDM e RSDM = —1.39. Later on, Amore|3] carried out a more 
accurate calculation with 3600 hydrogen eigenfunctions and obtained e RR = —0.128 
as well as e RRS = —0.132 by fitting and extrapolating the outcome of a Shanks 
transformation. The authors of both articles resorted to a variational parameter 
(decaying parameter [T| or length scale [3]) that considerably improves the result. The 
remaining disagreement between RR and RSDM is probably due to the well known fact 
that the basis set used in those calculations is not complete because it does not include 
the continuous spectrum[2] (see, for example Ref. jl] and the references therein). For 
this reason the RR calculation proposed by those authors [U |3] is expected to have a 
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limited accuracy no matter how large the dimension of the basis set of discrete states. 
It is surprising, however, that the lack of the continuous wavefunctions appears to be 
more noticeable for the ground state. In addition to it, at first sight it seems that there 
is a better agreement for the odd states^. 

The purpose of this paper is to carry out a RR calculation with a nonorthogonal 
basis set of square-integrable functions that in principle does not require the continuous 
spectrum. In section [2] we describe the RR variational method with such an improved 
basis set. In section [3] we compare present results with those obtained earlier by 
Dasbiswas et alpQ and AmorejH]. Finally, in section H] we summarize the main results 
and draw conclusions. 

2. Rayleigh-Ritz variational method 

The linearized model for the Ginzburg-Landau theory leads to the Schrodinger equation 

ti 2 ^9 , cos 9 „ , 

- 7^— V 2 ^ + P i> = Eip 1 

2m r 

where r = yfx 2 + y' 2 , < 9 = arctanQy/a;) < 2ir, and p is the strength of the dipole 

potentialpQ. Choosing the units of length h 2 /{2mp) and energy 2mp 2 /h 2 we obtain the 

dimensionless eigenvalue equation 

^9 , cos# , 
r 

where e = h 2 E / \2mp 2 ) . Since the potential V(r,9) = cos9/r is invariant under 
reflection about the x axis V(r, —9) = V(r,9), then the wavefunctions are either even 
ip(r, —9) = if;(r, 9) or odd ip(r, —9) = —ip(r, 9) under such coordinate transformation. 

Dasbiswas et alp] and Amore[3] resorted to the RR variational method with the 
basis set of eigenfunctions of the planar hydrogen atom: 
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where n = 1,2,... and R n> i(r) is the normalized solution to the radial equation[2]. 
However this basis set is incomplete if one does not include the eigenf unctions for the 
continuous spectrum (see, for example, [1] and references therein). 

If, on the other hand, the chosen basis set of square-integrable functions does not 
require the continuous spectrum to be complete then one expects the accuracy of the 
RR variational results to be determined only by the basis dimension. Here we propose 
the nonorthogonal set of functions 

j = 1,2,...} = {e~ ar , r l+1 cos' 9e~ ar , i, j = 0, 1, . . .} (4) 

for the even states and 

{4>°, j = 1,2, . . .} = p +1 sinfl cos'' 9e~ ar , i, j = 0, 1, . . .} (5) 

for the odd ones, where a > is a variational parameter. This basis set resembles 
the Slater orbitals commonly used in quantum chemistry calculations of atomic and 
molecular electronic structure [1]. 

The RR method with the variational ansatz 

j=m 

leads to the generalized eigenvalue problem 

HC = eSC (7) 

where Hij = (4>i\H\(j)j) and Sij = 4>j)- Note that it is possible to obtain explicit 
expressions for both kinds of matrix elements in terms of the variational parameter a; 
even more important is the fact that, with the help of a computer algebra software, like 
Mathematica, we can obtain the inverse of S explicitly for all the cases considered in 
the present paper. In this way the inversion does not introduce any round-off errors 
and the original generalized eigenvalue problem is converted to the ordinary eigenvalue 
problem 

S^HC = eC (8) 
for the nonsymmetric matrix S _1 H. 
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3. Results 
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For brevity we write both the exact and approximate RR eigenvalues and eigenf unctions 
as ei < 62 < . . . and ipi, ip2, ■ ■ ■, respectively. We express the rate of convergence of the 
RR results in terms of the largest parameter K = i + j in the wavefunction expansion 
fl6]), where i and j are the exponents of r and cos# in either equation (j4j) or ([5]). Thus, 
for a given value of K there are N = (K 2 + K + 2)/2 basis functions of either even or 
odd symmetry. The optimal value of the nonlinear variational parameter a depends on 
both the chosen eigenvalue and the number of terms N in the wavefunction expansion 
06]). For example, Fig. [TJ shows that for the ground state a increases with K oscillating 
about a straight line. Fig. |2] shows the RR eigenvalue e\ for a range of values of K. The 
rate of convergence is remarkably greater than the one for the Coulomb basis set [3]. 

Table [I] shows the RR eigenvalues obtained with the nonorthogonal basis sets (J4]) 
and (jSJ) {^n B ) on the one side and with the basis sets of even and odd Coulomb functions 
© ( e n B ) on the other. We appreciate the following facts: the accuracy of present 
nonlinear basis set is always greater (e^ B < e^ B ) in spite of the fact that the RR 
calculations have been carried out with with N — 211 nonorthogonal basis functions 
(K = 20) and N = 3600 Coulomb functions[3]. The discrepancy is greater for the lowest 
states and, it is less noticeable for the odd ones. We are presently unable to provide 
a rigorous proof for the last two facts; however, we may conjecture that the omitted 
continuous spectrum is not so relevant in those cases where there is agreement between 
the NB and CB variational results. The number of digits in the entries of this table 
is dictated by comparison purposes and does not reflect the estimated accuracy of the 
calculation. 

Both the analytical calculation of the matrix elements and the analytical inversion 
of the matrix S are time consuming but we do them only once for all the states. On 
the other hand, the optimization of the variational parameter a for every state is a 
time consuming calculation that we should repeat several times. Earlier and present 
calculations suggest that the RR variational method is less efficient for the ground state. 
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For this reason we have attempted variational calculations with considerably greater 
basis sets only for this state. For example, we have obtained e^ B even = —0.13774677227 
with N = 466 NB functions (K = 30) and e NBeven = -0.13774778205 with N = 821 
ones (K = 40). These results suggest that the first 6 digits remain stable and, 
consequently, that the RR calculation with the Slater-type basis functions does not 
appear to approach the RSDM results any more closely. However, it is worth noticing 
that present RR eigenvalues agree with the RSDM ones within the 2% error estimated 
by Dasbiswas et alp]. 

By means of the approximate ground-state wavefunction ipi(x,y) in terms of 
the Slater-like orbitals we have also calculated the effective dimensionless coupling 
constant [1] 

g = J dxdy^x^)] 4 (9) 

Dasbiswas et al[T] obtained g = 0.017 by means of a simple variational function 
constructed from the first three elements of the even nonorthogonal basis set (]1J): 
{e~ ar , re~ ar , r cos 8e~ ar } and g = 0.0194 by means of the RSDM. Amore[3] obtained 
g = 0.017 by means of a reduced RR basis set of Coulomb functions (N = 345). Note 
that just N = 3 NB functions yield the same result as N = 345 Coulomb ones. The RR 
method with N — 211 functions (K = 20) of the basis set (0} yields g = 0.0193 that is 
quite close to the RSDM result. 

Figures EJ |4] (left panels), [5] and |6] show the contour plots for t/)(x, y) 2 for the first 
five even and odd states obtained by means of the N = 211 NB functions. One clearly 
realizes that there are two types of nodal lines ip(x,y) 2 = and that the energy depends 
differently on each of them. The right panels of Figures |3] and H] show 3D plots of the 
probability densities of the first even and odd states. 

4. Conclusions 

Present results clearly show that the basis set of Slater-type orbitals is preferable to the 
Coulomb basis set. With just a few functions of the former set one obtains results that 
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are considerably more accurate than those arising from much larger sets of the latter. As 
argued above, the reason is that any linear combination of discrete-spectrum Coulomb 
eigenfunctions is orthogonal to the continuous-spectrum eigenf unctions. On the other 
hand, no continuous-spectrum functions are required when using the Slater-type basis 
set. The contribution of the continuous spectrum appears to be more relevant for the 
lowest states and for the even ones. However, we have not proved that the RR variational 
method with the Slater-type functions converges towards the actual eigenvalues as K 
increases. Present variational eigenvalues converge to limits that are slightly larger than 
the RSDM ones and we cannot safely state that all the stable digits of our results agree 
with those of the exact eigenvalues. Accurate lower bounds are required for that purpose 
and we have not yet been able to obtain them. In spite of this fact, it is encouraging 
that present RR results agree with the RSDM ones within the reported 2% accuracy 
of the latter pQ. If, as argued by Dasbiswas et alpQ, the RR variational method is 
more convenient than the RSDM for highly excited states, then present contribution is 
relevant because there is no doubt that the basis set proposed in this paper is preferable 
to the Coulomb one. 
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Table 1. Optimal a and energies for the first 5 even and odd states for K 
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Figure 1. Optimal value of a for the ground state as a function of K 
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Figure 2. Convergence of t\ in terms of K for the optimal value of 
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Figure 3. Contour and 3D plot for l^il 2 
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Figure 4. Contour plots for the squares of the second, third, fourth and fifth even 
states 
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Figure 5. Contour and 3D plot for the square of the first odd state 




Figure 6. Contour plots for the squares of the second, third, fourth and fifth odd 
states 



